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The existence of maximal and minimal mixing angles in the neutrino mixing matrix motivates 
the search for extensions to the Standard Model that may explain these angles. A previous study 
(C.I. Low and R.R.Volkas, Phys.Rev.D68,033007(2003)), began a systematic search to find the min- 
imal extension to the Standard Model that explains these mixing angles. It was found that in 
the minimal extensions to the Standard Model which allow neutrino oscillations, discrete unbroken 
lepton family symmetries only generate neutrino mixing matrices that are ruled out by experiment. 
This paper continues the search by investigating all models with two or more Higgs doublets, and an 
Abelian family symmetry. It is found that discrete Abelian family symmetries permit, but cannot 
■ explain, maximal atmospheric mixing, however these models can ensure #13 = 0. 
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Q H | zero in Eq. QJ. A special case of Eq. (JTJ is tri-bimaximal mixing, when sin 9 12 = 4j EJEIEIEJj and 6*13 is 
exactly zero. Two out of the three angles in Eq. Q assume extreme positions in parameter space - the minimum 
possible value and the maximum possible value - so it has been suggested by many authors that this mixing pattern 
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I. INTRODUCTION 



The approximate form of the neutrino mixing matrix 
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has been determined by neutrino oscillation experiments (Majorana phases have not been included). 

The mixing matrix, parameterised in the usual way, is formed by three very different mixing angles. The atmospheric 
mixing angle 623 has a maximal value of 7r/4 at best fit the solar mixing angle 612 has been found to be 

large 812 ~ 33° || , but not maximal, by solar neutrino oscillation experiments 0, 0, 0, El El , and the angle 
13, measured by the non-observation of v e disappearance |l3j |. is small and has only an upper bound and is set to 



is not accidental, but could be due to a family symmetry. 



' A. Family symmetry models 



The symmetries of the Standard Model (SM) do not dictate the Yukawa coupling strength between each fermion 
and the Higgs field. As a result, in the SM the charged lepton and neutrino mass matrices are 3x3 matrices with 
each element a free variable. In the SM, Dirac mass matrices have nine free variables, and Majorana mass matrices 
have six. Diagonalising the mass matrices generates the mixing matrix which can be any unitary 3x3 matrix. Family 
symmetries constrain the form of the neutrino and charged lepton mass matrices by relating elements of the mass 
matrix, or forcing elements to be zero, thus reducing the number of free variables in the mass matrix. 

For a family symmetry to fully predict a mixing matrix, all mixing angles must be independent of the free variables 
in the mass matrix, and must be prescribed by the form of the mass matrix. Mass matrices that can generate mixing 
matrices in this approach have been called "form-diagonalisable" matrices [Tsj . This can happen when there are three 
variables in the mass matrix corresponding to three unknown masses, and no free variables remaining for the mixing 
angles. An example of a form-diagonalisable matrix is the circulant matrix which can be generated by a Z3 symmetry 
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is diagonalised by 
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where to = e l27T ^ 3 and has eigenvalues a + b + c, a + tub + uj*c, a + co*b + ujc. 

The mixing matrix (Eq. |T)l) may be created by partially form-diagonalisable matrices, where the zero and maximal 
mixing angles are not related to any free parameter, and arise from the form of the mass matrices, while the solar 
mixing angle may be related to a free parameter. 

An Abclian family symmetry - individual lepton number U(1)l c x x U{\)l t - is conserved when neutrinos 

are massless, but is broken when the three neutrinos gain different mass values. The special form of l|T]l. and even 
more so the particular tri-bimaximal case makes it conceivable that a remnant of this Abelian group remains unbroken 
with massive neutrinos, constraining the mixing pattern. This further motivates the study of family symmetries, and 
Abelian symmetries in particular. 

Many models with family symmetries have been proposed. A number of these models |3 El produce the 
desired form of the mixing matrix, but use symmetries that cannot be easily incorporated into the SM as the left- 
handed neutrinos transform in a different way to the left-handed charged leptons, thus breaking SU{2)l- Most models 
that do preserve SU(2)l require additional fields such as singlet or triplet Higgs fields pH I22I V2^ or additional heavy 
fermions [24| • The models with the least new particle content require a number of Higgs doublets [2^, l26l |27| , and 
some soft symmetry breaking terms to generate Eq. l|T)l. 

It is clear that there are models that can produce Eq. (JIJ, the question that this work addresses is what is the 
minimal predictive model. The approach taken is to construct the simplest model, and find out whether the model can 
generate the mixin g m atrix, or whether it can be ruled out. If it is ruled out, the next simplest model is investigated. 
A previous study [18( began a systematic search to find the minimal extension to the SM that can generate the 
mixing matrix form of Eq. and found that for models with one SM Higgs doublet unbroken discrete Abelian 
family symmetries cannot produce the matrix. In fact these symmetries can only generate mixing matrices that are 
ruled out by experiment, or mixing matrices that are completely unconstrained by the symmetry. Non- Abelian family 
symmetries are also ruled out as they dictate that the charged leptons are degenerate. The structure of the next 
simplest model is a subjective question. I chose to study extensions to the SM with two or more Higgs doublets that 
transform under an Abelian family symmetry. Abelian symmetries were chosen as the symmetry group could be a 
subgroup of U(l)i, e x U(1)l x U(1)l t , and for simplicity. This case also differs from the single Higgs doublet case 
as it is possible for the exact family symmetry to be spontaneously broken by the Higgs VEVs. 



Section [H] presents the mathematics of family transformations, and shows how the mass and mixing matrices can 
be constrained by the family symmetry transformations. Three neutrino mass generation mechanisms are considered; 
left and right handed neutrinos coupling to the Higgs doublets to create a Dirac neutrino mass matrix, left handed 
neutrinos coupling to the Higgs to form a dimension-5 operator, and the seesaw mechanism where the right handed 
neutrinos get a bare Majorana mass. The types of mixing matrices that can be generated by an Abelian group are 
described in section ITTT1 I find that it is possible for Abclian symmetries to dictate that 6*13 — 0, and although the 
symmetries permit all atmospheric mixing angles, the symmetries cannot specify that the atmospheric mixing angle 
is maximal. Section llVI lists group transformation matrices that give a mixing matrix with #13 = 0. Section Ivl draws 
conclusions about these models and suggests other models and symmetries that may be more successful. 

II. HOW HIGGS AND LEPTON FAMILY SYMMETRIES CONSTRAIN MASS AND MIXING 

MATRICES 

The following section describes how a single transformation can restrict the Higgs-lepton coupling matrices. For 
a symmetry group of order n there are n of these transformations. However, if the Lagrangian is unchanged by a 



B. Outline 
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transformation X , it will also be unchanged by X m , where m is a positive integer. For Z n , the group of the addition 
of integers modulo n, the group is made up of the powers of one transformation, so a single transformation is sufficient 
to describe the restrictions placed on the coupling matrices. For all other groups more than one transformation is 
required. 



A. The symmetry transformations 



The family symmetry transformation matrices act on the different families of Higgs fields and leptons. The lepton 
transformation is 



Xr 



XtB.1} 



V R — * X uR V R , 



(4) 



where £l, vl, £r and v R are each 3-vectors in family space, £l and £ R are the vectors of left and right handed 
charged leptons, vl and v R are the vectors of left and right handed neutrinos. Each X matrix is a 3 x 3 unitary 
transformation matrix in lepton family space. To preserve SU(2)l, the left-handed neutrinos transform in the same 
way as the left-handed charged leptons. 
n families of Higgs fields transform via 



(5) 



where $ 



matrix. 



\<f>n J 



is an n-vector in Higgs family space, containing all the Higgs fields, and A$ is an n x n unitary 



B. Constraints on the Higgs- lepton coupling matrices from the symmetry 

1. The charged lepton Higgs coupling term 

The charged lepton-Higgs coupling transforms as 

h{^) T Xi R -^T L xl{^) T Al\X m £ Rl (6) 

where $° is an n-vector in Higgs family space containing just the neutral component of the Higgs doublet. The term 
(<1> ) T indicates a transpose in Higgs family space. A is an n vector in Higgs family space, where each element of A 
is a 3 x 3 Yukawa coupling matrix in lepton family space. Without any family symmetry the A matrices can be any 
3x3 matrices, but the existence of the family symmetry constrains them by 

A = x\Al\X m . (7) 

Note that Xt, commutes with $ as Xr, acts only on lepton family space, and $ is a lepton family singlet. The 
charged-lepton mass matrix is made up of the matrices in A: 

M e = (< $° >) T A. (8) 



2. The Dirac neutrino Higgs coupling term 
The Higgs field can couple to neutrinos in a number of ways. A Dirac neutrino coupling term transforms like 

VL($y KDiracVR ~» V : L Xl($ ) 1 A%K Dirac X vR V R , (9) 

where Krjirac is an n-vector in Higgs family space, made up of 3 x 3 Yukawa coupling matrices in lepton family space. 
Without the symmetry, the Yukawa coupling matrices can be any 3x3 complex matrices. Imposing the symmetry 
the matrices are constrained by 

^Dirac 

= X\A%K, mrac X vR , (10) 
which alters the Dirac neutrino mass matrix through 

M vDirac = (< $° >) ] n Dtrac . (11) 
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3. The dimension-5 neutrino Higgs coupling term 

A dimension-5 Higgs-neutrino coupling term transforms like 

jV l &k*'(v l )° - jV L Xl&Al K A%$*X* L (v L y, (12) 

where k is now annxn matrix in Higgs family space, with each component a 3 x 3 symmetric matrix in lepton family 
space, k is constrained by the symmetry through 

k = xlAl K A%X* L , (13) 

which consequently alters the mass matrix, defined by 

M„ = (< $° >)*«(< $° >)*. (14) 

4- Higgs coupling terms for seesaw neutrinos 

In the seesaw mechanism the right-handed neutrinos couple to form a bare mass term. The Higgs fields are not 
involved, so the mass term transforms as 

Vr^Mrvr -» Vr c XT r M r X vR v r . (15) 

restricting the heavy right-handed mass matrix Mr by 

Mr = X T vR M R X vR . (16) 

The resultant light neutrino mass matrix, given by M v — M Dirac M^ Mp irac is affected by the symmetry through 
the constraints on the heavy Majorana mass matrix and the Dirac neutrino mass matrix (listed in Sub. Sec. Ill B 2(1 . 
The seesaw case can be reduced to the dimension-5 operator case by relating 

KDira,cMRK^) irac — k from dimension-5 case. (17) 

^Dirac^lR^Di rac has all the constraints of k, plus additional restrictions from the transformation of the right-handed 
neutrinos. 

C. Abelian groups create mass matrices with zero or unconstrained elements 

The restrictions family symmetries have on mass matrices depend on the transformation matrices that are chosen. 
If the set of matrices Xi form a group then Yi = U^XiU form the same group, where U is any unitary matrix. Sets of 
matrices related in this way are called equivalent representations. Appendix^shows that choosing different equivalent 
representations for the lepton family symmetry transformations corresponds to choosing a different weak basis for 
the leptons. The constraints on the masses and on the mixing matrix are identical for two different equivalent 
representations. This makes it possible to eliminate groups, as each group has only a finite set of non-equivalent 
representations . 

This result simplifies the study of Abelian groups. As all Abelian groups are equivalent to a diagonal representation, 
only these representations need to be considered. Since the transformation matrices must be unitary, the diagonal 
elements are phases. 

This makes the charged lepton restriction of Eq. Q become 

Xij = {Xl)nA^Xij{X tR )jj no summation, (18) 

where i, j are lepton family indices. 

This restriction means that A}a (the ijth element of the Yukawa coupling matrix for <f>\) can be related to the 
ij element of the other Yukawa coupling matrices. The symmetry however, does not relate an element of A with a 
different element of A. If the ijth elements of all A matrices are zero then Muj = 0, otherwise Muj will most likely be 
non-zero, and unrestricted by the symmetry. However, if two or more ijth elements are related, and the Higgs VEVs 
are related then there could be cancellation: M^j — Kj < ( I >0 >= 0. A relationship between VEVs is possible as the 
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symmetry also constrains the form of the Higgs potential. If there isn't a cancellation the element of the mass matrix 
is unrestricted by the symmetry - it is a free parameter of the model. 

Consequently, a symmetry does not dictate the relationship between any elements in the mass matrix. What the 
symmetry does do is force some elements of the mass matrix to be zero, leaving all other elements unrestricted. This is 
true for instances where there is cancellation, and when there isn't a cancellation. This makes the analysis of Abelian 
groups easier, as only mass matrices with zero and unrestricted elements need to be considered, and analysis of the 
Higgs potential is not required. 

There is only one Dirac mass matrix of this type that is form-diagonalisable - the diagonal matrix which is 
diagonalised by the identity - so most charged lepton diagonalisation matrices will depend on the elements that are 
unrestricted by the symmetry. Partially form-diagonalisable matrices are possible, for example a mass matrix which 
is in 2 x 2 block diagonal form, is diagonalised by a unitary matrix which is in 2 x 2 block diagonal form, which has 
one Euler angle depending on the free parameters, and the other two angles zero. 

Similarly, the Dirac neutrino-Higgs coupling matrices are restricted by a diagonal transformation, 



K-Di 



(xt: 



i-^^^Dirac ij {.-^vfljjj : 



(19) 



which also yields mass matrices with elements that are either zero or unrestricted. The dimension-5 neutrino-Higgs 
coupling matrix is constrained by 



(X^ L )uAl,K i jA^,(X^)j 



(20) 



which also constrains some mass matrix elements to be zero. However, as Majorana mass matrices are symmetric, 
more form-diagonalisable mass matrices can be created. These mass matrices have a pseudo-Dirac form: 




is diagonalised by 
is diagonalised by 
is diagonalised by 
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(21) 
(22) 
(23) 



The neutrinos that arc mixed have 



Partially form-diagonalisable matrices also can be created 
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(24) 



(25) 



(26) 



where 9 and a are angles which depend on the parameters a and b. Again, the two neutrinos that are maximally 
mixed have m, = m, . 



Right handed Majorana mass matrices are constrained by a diagonal transformation by 



(27) 



and can also generate the form-diagonalisable matrices in equations l(2*Tl - l2*6" |l . 
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D. Non-Abelian groups 

For family symmetries where the Higgs fields do not transform, non-Abelian symmetries ensure that at least two 
charged leptons must be degenerate 18]. In this case the mass matrix is constrained by the symmetry through the 
equation Mi — X^MiXir. Consider the mass basis, where Mi is diagonal. Non-Abelian representations cannot be 
equivalent to a diagonal representation, so the transformation associated with the mass basis will mix the mass matrix 
elements, and ensure that at least two of the masses are equal. For cases with more than one transforming Higgs field, 
the mass matrices are made up of a number of Yukawa matrices, and the transformations act on these matrices in a 
more complicated way than the one Higgs field case. Non-Abelian transformations no longer neccesarily force equal 
mass constraints, and as a result cannot be ruled out in the same way as in Ref |l8j . 

III. MIXING ANGLES THAT CAN BE GENERATED BY ABELIAN GROUPS 

A. Mass matrices to investigate 

To find the types of mixing matrices that can be created by Abelian groups a program was created to generate all 
sets of neutrino and charged lepton mass matrices with zero and unrestricted elements. For a given set of neutrino 
and charged lepton mass matrix types two sets of mass matrices were created - each with the same textures (i.e. the 
same positions of the zeros), but different random numbers were used for the elements that were unrestricted. The 
unitary diagonalisation matrix was found for each mass matrix, and two mixing matrices were created and compared. 
If an angle was the same for both mixing matrices, the value of the angle was due to the textures of the mass matrices, 
and thus, a result of the symmetry. If an angle was different for the two mixing matrices then the angle's value was 
due to the random numbers in the mass matrices, and not concerned with the symmetry 

This was done for Dirac neutrinos, Majorana neutrinos which gained mass from a dimension-5 operator, and seesaw 
neutrinos. For seesaw neutrinos only right-handed Majorana neutrino mass matrices that were invertible were used, 
as non-invertible matrices generate less than three ultra-light neutrinos [28l 12^ . 

B. Results 

1. Abelian groups can generate zero and maximal mixing angles 

The only form-diagonalisable Dirac mass matrices that can be generated by Abelian groups create diagonalisation 
matrices with Euler angles equalling zero (Sec. Ill C() . Consequently, mixing angles for Dirac neutrinos are either zero, 
or unfixed by the symmetry, meaning Dirac neutrino models can ensure #13 = but cannot fix the atmospheric or 
solar mixing angles. These angles will be free parameters and can take any value. In some cases these two mixing 
angles can be related to lepton masses. 

Majorana neutrino mass matrices, from seesaw and dimension-5 operators, can create fixed zero mixing angles, 
and also create fixed maximal mixing angles from the pseudo-Dirac type mass matrix (Eq. Q21l - I21)|) ). This looks 
promising for creating a maximal atmospheric mixing angle, and the fact that maximal mixing angles can only be 
generated from Majorana matrices perhaps could be a key to explaining why lepton mixing angles are large but quark 
mixing angles are not. Unfortunately, it was found that the maximal mixing angle cannot correspond to atmospheric 
mixing. 

2. Conditions for maximal mixing 

Mixing angles that are fixed to be maximal can only arise from the form-diagonalisable Majorana mass matrices 
listed in equations 1211 - 1261 These matrices have a 2 x 2 pseudo-Dirac block and the maximally mixed neutrinos 
always have m, = —rrij, corresponding to Sm 2 = 0. This means that there will be no oscillation, however a small 
mass squared difference could be created by breaking the symmetry. 

To demonstrate the difficulty in generating maximal atmospheric mixing, consider mixing matrices that have #13 = 
and one maximal mixing angle. Matrices of this type must have a pseudo-Dirac neutrino diagonalisation matrix (Eq. 
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Q21I— I26[l) and a charged lepton diagonalisation matrix that has either zero or one mixing angle. That is, 
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where a and <5 are angles undefined by the symmetry. 

Charged lepton diagonalisation matrices with more than one mixing angle cannot produce mixing matrices with 
zero and maximal mixing angles, so the only combinations of diagonalisation matrices that have both #13 = and a 
maximal mixing angle are 
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where the U„i's are given in Eqs: 1)21126(1 

The only mixing matrix that has maximal atmospheric mixing is E/mnS2, which also has a very unsatisfactory solar 
mixing angle of zero. C/mnS3 has maximal solar mixing, but the atmospheric angle is not dictated by the symmetry. 

The program that was written (see Sub. Sec. IIII Afl also searched for fixed mixing angles without the 0i3 = 
constraint and found that # 2 3 is still unfixed by the symmetry, therefore the only aspect of the mixing matrix form of 
Eq. UJ that can be generated by a symmetry is #13 = 0. It is not possible to demonstrate this result in a concise way 
in this paper: it is instead the result of a systematic computer-aided search. We have seen that the mixing angles 
can be zero, maximal or unfixed by the symmetry. The unfixed mixing angles can be either related to fermion masses 
or completely free variables. The possible ways in which mixing angles can be related to masses was not analysed 
by the program, however, for the cases where #13 = mass-mixing angle relationships were worked out by hand 
rSub.Sec HTBl . 

Note that the fact that solar mixing can be forced to be maximal and atmospheric mixing angle cannot be forced to 
be maximal does not indicate any fundamental difference between the flavours. If one mixing matrix can be predicted 
by an Abelian group, so can the mixing matrix with rows permuted. Permuting the rows corresponds to interchanging 
e, fj, or r, so the whole set of possible neutrino mixing matrices are flavour symmetric. 

The mixing angles, defined as Euler angles, however are not flavour symmetric. The probability that a neutrino of 
flavour i is detected as flavour I' after a distance x is given by 
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(32) 



where L mm i — 2tt ■ When #13 = 0, the probability of an electron neutrino being detected as an electron neutrino 

after a distance x is 
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and is only dependent on 0i2. The probability of a muon neutrino being detected as a muon neutrino after a distance 
x, 
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(34) 



(35) 



is dependent on both non-zero mixing angles. So when 0i3 = a maximal solar mixing angle corresponds to a 
maximum amplitude of oscillation - an electron neutrino will oscillate into a state with no electron neutrino component. 
Maximal atmospheric mixing means that there is a mass eigenstate that is an equal superposition of Vn and v r and 
does not imply a maximum amplitude of oscillation in the three flavour case. 
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TABLE I: Mass matrices for Majorana neutrinos that give #13 = 0. Two Higgs doublets are required unless otherwise stated. 

IV. MASS MATRICES AND SYMMETRIES THAT PRODUCE 6 13 = 

There are several sets of charged lepton and neutrino mass matrices that can produce #13 = 0, some of which 
can be created from a symmetry. The ones that can be related to a symmetry, and do not force the muon or tau 
leptons to be massless are listed in tables III and ITT1 along with the smallest symmetry group that can produce the mass 
matrices. All cases that can be generated by a symmetry require two Higgs doublets, unless otherwise stated in the 
table. Cancellation within the mass matrix was not considered, nor was the possibility of VEVs equalling zero (i.e. 
if Mg 3 = 0, then it was assumed that \ % i 2 n = 0)- With these assumptions, diagonal Higgs transformations give 
the same mixing matrices as equivalent non-diagonal representations (see App. 0), so to find the smallest symmetry 
group only diagonal transformations were investigated. It is possible that a smaller group than that listed could 
produce the mixing matrices if there are cancellations or zero VEVs. 

The smallest group that can give 6*13 = is Z3, the group of addition modulo 3. Z2 gives the same mixing 
matrices that can be generated in the single Higgs doublet case, as analysed in Ref. |1J| and are either unrestricted 
or experimentally ruled out. This is shown in App. [C] 
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TABLE II: Mass matrices for Dirac neutrinos that give #13 = 0. Transpositions of the columns of the mass matrices do not 
alter the masses, the mixing matrix, or the symmetry. Two Higgs doublets are required unless otherwise stated. 



A. Examples of symmetry transformations 

For seesaw neutrinos the mass matrices in the first row of table [I] can be generated by the Z 4 transformation 
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The first row of mass matrices in table |n] can be generated by the Z 4 transformation 




(36) 
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(37) 



For many cases the Dimension-5 mass matrices require larger symmetries or additional Higgs doublets. This is 
because the seesaw and Dirac cases have extra freedom due to the transformation of the right-handed neutrino. 
The first row of matrices in table [i] can be generated by a Z-j transformation involving three Higgs doublets: 
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where (3 = e 2lxl l~ ! . 
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B. Mass and mixing angle relationships 



Section IIII Bl established that none of the mass matrices that give #13 = can give a fixed atmospheric mixing 
angle. In fact, it can be shown by diagonalising the mass matrices of tables [I] and |nj that the atmospheric mixing 
angle is also unrelated to masses - it is a free variable for all #13 = cases. 

In cases where a Majorana mass matrix is the source for the solar mixing angle, the mixing angle can relate to the 
masses. These cases are the Majorana neutrino matrices from rows 2-4 and 5 - 8 of table [I] 

The neutrino mass matrix of rows 2 and 6 of table [I] 



relate the masses and mixing angle by 



/ 


.4 


B 







B 








V 








D 



(39) 



tan 26 



2 A /-mim 2 
m i + m 2 



(40) 



Using the approximate values for 8m\ 2 ~ 7.5 x 10 eV and so i w 33°, the neutrino masses must be approximately 
hierarchical, with \mi\ ~ 9 x 10~ 3 eU, |m 2 | ~4x 10~ 3 eU, and |m 3 | ~ 0.04 eV. Rows 4 and 7 give the same mass 
pattern except mi and m 2 are interchanged. 

The mass matrices of rows 5 and 8 are even more constrained. There are only 3 free variables in these neutrino mass 
matrices, and these free variables describe two masses, the solar mixing angle, and a contribution to the atmospheric 
mixing angle. As a result the mixing angles and masses must be related: 



mi 

m 2 
m 3 

tan(20i 2 ) 
tan 023u 



i(| J 4|±V|A| 2 +4(|B|2 + |C|2 ))! 



-(\A\ T vW+WF+m), 



= 0, 



A 



\C\ 
\BV 



(41) 

(42) 
(43) 

(44) 
(45) 



where (9 2 3„ , along with the diagonalisation angle from the charged lepton mixing matrix form the atmospheric mixing 
angle. To achieve 5m\ 2 "C Sm 23 , \A\ <C -\/|-E>| 2 + |C| 2 is required, and the neutrino mass matrix becomes close to the 
partially form-diagonalisable matrix of Eq. 1241 giving very nearly maximal solar mixing. This has been ruled out by 
experiment. 



C. Flavour changing neutral currents 



Models involving a number of Higgs doublets and Abelian symmetries naturally predict flavour changing neutral 
currents (FCNCs) for charged leptons and neutrinos j2jj, and the models presented in this paper are likely to be no 
exception. However it is possible that the flavour symmetry somehow suppresses the FCNCs. Many of the charged 
lepton mass matrices listed in tables HI and ITT1 only mix /i and r, and do not mix electrons, meaning that the most 
experimentally constrained FCNC processes (such as /i — ► eee) are not allowed. However the decay r — > Ji/J,/! is 
allowed, and as large off diagonal elements in the Yukawa coupling matrices are required to give large mixing angles, 
this transition is not likely suppressed. However, it is possible that some action of the flavour symmetry on the Higgs 
fields can prevent the FCNCs from becoming too large. 



V. CONCLUSION 



The best fit neutrino mixing matrix Eq. Q has a few peculiar aspects; it is very different from the CKM matrix, 
and it has one maximal mixing angle and one minimal mixing angle. It would be pleasing to find that this pattern 
can be generated by a symmetry. 
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Earlier work 18] showed that unbroken lepton family symmetries alone can only produce mixing matrices which 
are not allowed experimentally. This paper continues the search for a symmetry explanation to the form of the mixing 
matrix. The models considered are extensions of the SM that include a number of Higgs doublets and discrete Abelian 
symmetries that transform the Higgs and lepton families. 

Symmetries of this type can only fix mixing angles to be zero or maximal, otherwise the angle can be any value 
as it depends on the free parameters of the model. This is due to the fact that all Abelian representations are 
equivalent to diagonal representations. In the diagonal basis Abelian symmetries can only dictate whether an element 
in a mass matrix is zero or unrestricted; no relationships between mass matrix elements can be generated, so only a 
few form-diagonalisable mass matrices can be generated, and most mixing angles are not fixed by the symmetry. A 
small number of Majorana mass matrices can generate maximal mixing, however, this mixing cannot correspond to 
the atmospheric mixing angle. The characteristic of Eq. Q that can be produced by Abelian family symmetries is 
8 13 = 0. This requires at least two Higgs doublets and a Z% or larger family symmetry. 

Although Abelian symmetries have limited ability in predicting fixed mixing angles, symmetries can relate lepton 
masses and mixing angles. For the cases where #13 is forced to be zero, the solar mixing angle can be related to 
neutrino masses. This relationship fixes the neutrino mass pattern to be hierarchical. 

Non- Abelian family symmetries may produce better results, as they can relate different elements in the mass matrices 
together, possibly creating form-diagonalisable matrices that cannot be generated with Abelian groups. Extending 
the Higgs sector by including triplet Higgs fields to generate neutrino mass is also likely to increase the possible types 
of mixing matrices. 

The approach taken here, to find the minimal model that explains the neutrino mixing matrix, has succeeded 
in explaining one of the interesting aspects of the mixing. However, the models that can explain 6* 13 = are not 
particularly simple. Fixing this one variable requires the introduction of extra Higgs doublets, which additionally can 
lead to flavour changing neutral currents. These results suggest that the minimal model route may not readily yield 
a satisfactory explanation for the mixing parameters. This could be a consequence of considering neutrino mixing 
independently of other unresolved issues in particle physics, such as mass hierarchy and quark mixing. Perhaps the 
answer can only be found by finding a model that simultaneously addresses several of these problems. 
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APPENDIX A: EQUIVALENT REPRESENTATIONS OF THE LEPTON TRANSFORMATIONS YIELD 

IDENTICAL PREDICTIONS 

Two equivalent representations for the fermion transformations, Al, Am, A vR and Bl, Bg R , B vR , are related by 

B L = S\A L S L , (Al) 

Bm = S tR Ai R S^ R , (A2) 

B v r — Sl R A vR S vR . (A3) 
where Sl, Sg R and S vR can be any 3x3 unitary matrices. 



1. Charged leptons 



The Yukawa matrices arising from the B transformations are denoted by a B subscript. The restrictions from the 
B transformations on the charged lepton Yukawas (from Eq. (|1U|) ^ are 

\ B = B\Al\ B B m (A4) 

= SlA\S L Al\ B s\ R A tR S lR . (A5) 

Rearranging gives 

Sl\bS\ r = A\S L Al\ B S\ R A m (A6) 

= AlA\{S L \ B S\ R )A lR . (A7) 
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(Xj commutes with U^Sl, as they operate in different spaces). 

The charged lepton Yukawa restrictions for A transformations are A,4 = A^A^AaAir. SlA b SI r has the same 
restrictions from the symmetry as A a- As the mass matrices are completely unconstrained apart from the generation 
symmetry constraints, we can set SlAbSI r = Aa- This means that the charged lepton mass matrix from the A 
representation can be given by 

Mia = Aa < $ > T = S L A B S\ R < $ > T = S L M m S\ R . (A8) 
If MgA is diagonalised by Uf L , and Uf R , then Mg B will be diagonalised by Uf L = SjJU^ and Uf R = Sf R Uf R . 

2. Dimension-5 neutrino masses 

Similarly, the conditions on the dimension-5 neutrino coupling matrices mean we can identify 

ka = SlkbSI, (A9) 

and the two mass matrices can be related by 

M uA =< $ >t K A < ^ >*=< * > f S l k b SI < >*= S L < $ >t k b < $i >* Si = SlAUbSI. (A10) 

M v a is diagonalised by U vA , and M v b is diagonalised by U v b — SlU v b, giving identical mixing matrices for the 
two transformations, 

C/mnsb - Ugll* = U^S{S L U V = U MmA - (All) 
3. Dirac neutrinos 

The neutrino Yukawa restrictions for the two transformations are related by SL,KDiracBS vR — KDi rac A, and the two 
mass matrices can be related by SlM v bS}, r = M VA . This gives diagonalisation matrices related by U^ L = SlU^ l . 
The mixing matrix is, therefore, the same for both transformations: 

£^mns = V%U? L = U£& L S L U? L = U£ NS . (A12) 

4. Seesaw neutrinos 
The right-handed Majorana mass matrix is restricted by Eq. (|lfi|l . giving 

Mrb = B^ r M R bB vR = S^ r A£ r S* r Mr B S uR A v iiSur, (A13) 

while M R a is constrained by Mra = A^ r M R aA vR . M R a can be equated to S^rM R bS^ r , as they have the same 
constraints from the symmetry. 

The Dirac neutrino mass is as above: SLM B i rac BS^ R = Mr>iracA- The resultant neutrino mass matrix is 

M v a — M Dirac aM R x Mp iracA — S lM DiracB M R BM D . iracB S 1 [ = SlM vB S l . (A14) 

The seesaw mass matrices are related to each other in the same way as the dimension-5 mass matrices. Using the 
result from Sub. Sec IA 21 the mixing matrices from the two representations are equal. 

APPENDIX B: EQUIVALENT REPRESENTATIONS FOR THE HIGGS TRANSFORMATION GIVE 
THE SAME MIXING PREDICTIONS IN MOST CASES 



Two equivalent representations for the Higgs transformations A$, and B$, are related by 

-B$ = S^A^S. 



(Bl) 
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The Yukawa matrices arising from the A and B transformation are denoted by an A or B subscript. 
The charged lepton Yukawa restrictions for the B transformation are 

X B = X{BIX B X IR (B2) 

= x{S T AlS*X B X lR . (B3) 

Rearranging gives 

S*X B = xlAlS*X B X m . (B4) 



The restriction from the A transformation is 



xlAlX A X m . (B5) 



Xa has the same restrictions as S*Xb, so they can be equal: Xa — S*Xb- The mass matrix Mia is a linear 
combination of A^ 2 ' therefore it is also a linear combination of X g ,2, "' n . If the symmetries do not dictate the 

ratios between the Higgs VEVs (e.g. for two Higgs fields and are unfixed) then Mib can be any linear 

combination of As, and Mgs can be any linear combination of Xa- Therefore Mib has the same restrictions from the 
symmetry as Mia, and the different equivalent representations make the same predictions. 

If there is a relationship between the VEVs it is possible to get extra zeros in the mass matrices. This occurs when 
the ijth elements of the A matrices are non-zero, but MV =< 0° > T A lJ = 0, due to a special relationship between 
the VEVs and elements in the A matrices. Also if one of the VEVs from a particular representation is equal to zero, 
then it is also possible for more zeros to be created in the mass matrix. 

The neutrino mass matrices are also unchanged by a change of representation . Both M v a and M v b are linear 
combinations of the same k matrices. The dimension-5 Higgs-neutrino coupling has 

K B = x{bIk b BIXI (B6) 
= X[S^ A^SkbS 7 'A%S* X* l (B7) 
(Sk b S t ) = xIA${Sk b S t )AIX* Li (B8) 

so SkbS t can be equated to ka, and M v a and M v b are both linear combinations of Kb matrices. 
For Dirac neutrinos 

K B = x{bIk b X vR , (B9) 
= xl&JkSKBXuR, (BIO) 
(Sk b ) = xIa%{Sk b )X vR , (Bll) 

so Skb and ka can be equated, and M v a and M v b are linear combinations of the same n matrices. Again, if the 
ratios between the mass matrices are not defined by the symmetry, any linear combination is a valid mass matrix, 
therefore both mass matrices are constrained by the symmetry in the same way. 

When there is cancellations or zero VEVs, the only change to the mass matrix is some extra zero elements. 



APPENDIX C: Z 2 SYMMETRIES GENERATE MIXING MATRICES THAT ARE EITHER RULED OUT, 

OR UNCONSTRAINED 

Applying a Z2 Higgs transformation twice leaves the Higgs fields unchanged, so for the mixing matrix to be allowed 
by the no-go theorem of 0|, -^-1 = -^m = -^uR = ^ or diagonal Higgs transformations, the components of ^4$ will 
also be 1 or — 1 . This appendix shows that mass matrices generated by a Z 2 transformation have equivalent restrictions 
to mass matrices generated by family symmetries when the Higgs field is not transforming. These restrictions are 
Mg — X^ L MfXf R for charged leptons, Mu irac — x\MB> irac X vR for Dirac neutrinos, M u = x\M v X'^ for dimension-5 
neutrinos, and M R = X^ R M R X vR for right-handed Majorana neutrinos. These situations have been ruled out by the 
theorem in |l8j. 
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1. Assuming no cancellations and no zero VEVs 

Due to the result of appendix iBl diagonal Higgs transformations can be used, and can have 1 or —1 as diagonal 
elements. For A$ = I the Higgs fields do not transform - this is equivalent to the single Higgs field case. For A<$, = —I, 
the restrictions on the charged lepton Yukawas reduces from A = Xj^A^XXm to A = — X^XXtR - which is equivalent 
to a single Higgs doublet scenario where the right-handed Higgs fields transform with —X^r. The dimension-5 neutrino 
restrictions do not change and the Dirac neutrino Yukawa restrictions are the same as the single Higgs case when the 
right handed neutrinos transforms under — X v r. 

When A® is made up of both Is and —Is, the restrictions on the charged lepton mass matrix are (from Eg. 1118(1 ) 
Mf = unless (x\) ll Xf R = 1 or -1. 

This condition holds for all i and j, so the charged lepton mass matrix is unrestricted by the symmetry. The 
restrictions on a dimension-5 neutrino mass matrix are similar (from Eq. 1)20(1): M l j = unless = 1 or 

— 1, which also will hold for all i and j, meaning that the neutrino mass matrix is also unrestricted by the symmetry. 

Dirac neutrinos will also be unrestricted, and as a result, the seesaw neutrinos will be unrestricted by the symmetry. 

2. Including the possibility of cancellations 

a. Charged leptons 

Cancellations in the charged lepton mass matrix means that for some i, j, M^j — A y - < $i > +Xfj < $2 > +■•• = 0, 
while the Xij are non-zero. This will only occur for particular ij's where there is a certain relationship between the 
A's, and the VEVs. For Z2 Higgs transformations there can be only two possible relationships: 

• Ay = A^Xij, which occurs when {X\) U X^ R = 1. If there is a cancellation for this Yukawa relationship, the 
cancellations add zeros into the mass matrix in exactly the same way as the condition Mi = —X* L M^XtR. 

• Ay- = — A^Xij, which occurs when = — 1. If this relationship led to a cancellation, the restrictions 
on the mass matrix would be identical to the conditions from Mi = X' L MtXiR. 

b. Dimension-5 neutrinos 

If there is cancellation in the neutrino mass matrix, the cancellation will occur for one of two relationships between 
the coupling matrices. 

• Cancellation when = A^k^A^, is equivalent to the restriction M v = — X\M V X*^. 

• Cancellation when k 1 -? = — A\k^ A%, is equivalent to the restriction M v = X\M V X^- 

All combinations of charged lepton mass matrix restrictions and neutrino mass matrix restrictions are the same as 
restrictions for single Higgs field cases. 

c. Dirac neutrinos 

The two relationships between the Yukawa coupling terms are — ztA^K^ . 

• Cancellation when k 4 - 7 = is equivalent to M v = —X^ l M v X v r. 

• Cancellation when k 4 - 7 = — is equivalent to the restriction M v = x' i L M v X V R. 

All combinations of neutrino and charged lepton mass matrix restrictions are equivalent to single Higgs cases, in the 
same way as the dimension-5 neutrinos. 
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d. Seesaw neutrinos 

The cancellation in the Dirac neutrino mass matrix is just the same as above, the heavy Majorana mass matrix is 
restricted by the symmetry by Mm = XI r MmX* r . The constraints from all cancellation possibilities arc identical 
to restrictions from single Higgs doublet cases . 

3. Zero VEVs 

If the Higgs transformation is diagonal, and < (j>\ >= 0, then the mass matrix is a linear combination of A 2 ' 3 '"'™, 
which obey A^X^X 01 X^r — X a , where a is the Higgs family index ranging from 2 to n. These restricitons are identical 
to a case with one fewer higgs field. 

For a non-diagonal Higgs transformation, the different A matrices are related by A = X^A^XXgR. For Z 2 transfor- 
mations this reduces to Ay = A J Ay if (XjJaXfy = +1, and Ay = —A J Ay if (X^uXfy = — 1 - these are just two 
sets of simultaneous equations. If the equations allow one or more of A 2 ' 3 --" to be non-zero, the M^j is non-zero, 
otherwise the element of the mass matrix will be zero. 

There arc 4 possibilities: 

• Muj = if (X\)nXtRjj = 1 otherwise Muj is unrestricted. This is the same restriction as Mg = —X\MiX^r. 

• Muj = if (X^uXiBjj = —1 otherwise Muj is unrestricted. This is the same restriction as Mg — -\-X\MgXgR. 

• M ti j = for all (Xl) u X mjj . 

• Maj is unrestricted for all (X' L )aXgRjj. 

For Dirac neutrinos the situation is similar. The KDirac matrices are related by KDirac = X* l A'^XX v r. The same 
four possible restrictions arise: Moirac = J rX\ i Mr>i rac X V R, Moirac = —X^MmracX^R, Mmrac is unrestricted, and 

Moirac = 0. 

For neutrinos with masses due to a dimension-5 operator, the n matrices are related by k = X^A^kA^X^- The 

possible restrictions on the neutrino mass matrix are: M v — +x|M„X£, M„ = +x\ j M v X'L M v is unrestricted, and 
M v = 0. 

A right-handed Majorana mass matrix is unaffected by the VEVs, so the usual restriction applies: Mr = 

X„rMrX„r. 

The combinations of neutrino and charged lepton mass matrix restrictions give four possibilities: The restrictions 
are identical to single Higgs field cases, the neutrinos are massless, the charged leptons are massless, or the mixing 
matrix is unconstrained by the symmetry. Therefore a Z 2 transformation predicts mixing matrices that are either 
ruled out or unconstrained by the symmetry. 

Note that this does not mean that groups which have Z2 as a subgroup can be ruled out. If the Z2 subgroup gives 
mixing matrices that are not allowed, then the group is ruled out. However, if the Z 2 subgroup leaves the masses and 
mixing angles unrestricted, then the group is still allowed. 
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